1 



as 

Oh! 
D 

(N 

> 
00 

O 
0^ 
O 



Oh! 



X 



RU-94-72 

Large Gauge Theory - Expansions and Transitions 
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We use solvable two-dimensional gauge theories to illustrate the issues in relating large A'^ gauge theory to 
string theory. We also give an introduction to recent mathematical work which allows constructing master fields 
for higher dimensional large A'^ theories. We illustrate this with a new derivation of the Hopf equation governing 
the evolution of the spectral density in matrix quantum mechanics. 

Based on lectures given at the 1994 Trieste Spring School on String Theory, Gauge Theory and Quantum 
Gravity. 



1. Introduction 

The idea that QCD can be reformulated as 
some sort of string theory refuses to die. If we 
allow ourselves a sufficiently liberal definition of 
'string theory' - I will take this to be 'a theory 
of embeddings of two-dimensional manifolds into 
space-time with an action local both on the world- 
sheet and in space-time,' allowing other degrees 
of freedom as well and noting that the possible 
theories have not yet been classified - we should 
agree that the idea has not been disproven. Nev- 
ertheless it has been difficult to make progress 
with it. 

In these lectures I will describe some of the 
work which has been done in the last two years 
on two dimensional gauge theory. Although this 
is drastically simpler than the four dimensional 
theory and can only illustrate a few qualitative 
features common to both cases, some important 
ideas have come out of this work, which suggest 
new lines to pursue in four dimensions. 

A good starting point for discussing 'QCD 
string' is the strong coupling expansion, which 
I will review very briefly. (See the large N 
case was discussed in detail by V. Kazakov and I. 
Rostov at the 1993 Trieste Spring School |3|,|3|.) 



We write the path integral 
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expand in a power series in /? 1/5^, integrate 
termwise, and sum. To do this, it is necessary to 
be able to make sense of the functional measure 
DA without benefit of the Boltzmann weight, and 
so far the only way known to do this in D > 2 is 
to latticize the theory. We thus use the Wilson 
action S[U] = J2p U{P)+U+{P) and link inte- 
grals such as JdUiL) U{L)lU{L)+l = 1^5^^5^^ 
to evaluate the 0(/3") term as a sum over dia- 
grams with n plaquettes glued together in a con- 
tinuous and 'surface-like' way. This expansion 
has a finite radius of convergence (because the 
number of diagrams only grows exponentially in 
n) and clearly produces area law behavior for Wil- 
son loops at small (3. 

At large N, one can interpret the rules for 
building diagrams in a way which makes a dia- 
gram with weight N'^~'^9 into a continuous two- 
dimensional manifold of genus g. Basically, this 
is because each link integral, contributing an edge 
to the diagram, produces 1/iV; the usual large N 
coupling dependence gives each plaquette an A^, 
while each vertex comes with a sum over an in- 
dependent index, giving N . (The roles of vertices 
and faces are switched from the weak coupling ex- 
pansion.) The main subtlety is dealing with link 
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Figure 1. A world-sheet 'saddle' associated with 
terms like (3). Every edge is embedded in the 
same lattice link. 



integrals such as 

J ]l 32 33 3i 33 31 3i 32 ^ ' 

33° 32° 3i" 31 ^''> 
+3 ^ 4. 

The subleading terms in this expression can con- 
tribute to a diagram with overall weight N"^ , and 
thus we must interpret an appearance of this term 
as a geometric object which can be glued to the 
rest of the surface to form a sphere. This will 
work if we call it a 'disk' with four edges and zero 
area, as in figure 1 ~ the geometrical TV count- 
ing is 7V^~'' which agrees with the result from 
the integral. The minus sign is to be interpreted 
as a weight associated with this feature. Higher 
'saddles' with more edges also appear in the full 
expansion. 

This already looks pretty close to the defini- 
tion of 'string theory' I started with and thus the 
project of finding a world-sheet action associated 
with some continuum limit of this to get a 'gauge 
string' seems very well motivated. Already on 
this qualitative level, many differences with fun- 
damental string theory are apparent. Perhaps the 
most striking is that the world-sheets are contin- 
uously embedded in the target space: this is in- 
compatible with a world-sheet action / d?(j{dxy' 
and not easy to describe with a unitary world- 



sheet theory, but on the other hand it may explain 
why in QCD string theory (unlike the fundamen- 
tal string) we can define correlation functions of 
operators local in target space. 

A serious objection can be raised to start- 
ing from a strong coupling expansion: to take 
the continuum limit of the theory, we must take 
(7^0 and the lattice spacing to zero as pre- 
scribed by the RG, and surely it is unrealistic to 
hope that the expansion will have infinite radius 
of convergence. One response to this objection 
is to argue that our string theory will agree with 
gauge theory only at long distances, and that this 
is an acceptable limitation. I disagree with this 
response. On theoretical grounds, if this is what 
we want, we can simply use the effective long- 
distance theory defined by Polchinski and Stro- 
minger p5[ , because no structure specific to the 
underlying gauge theory survives in the long dis- 
tance limit. On practical grounds, what is lacking 
is a precise calculational technique which works 
for intermediate scales, interpolating between the 
known short and long distance behaviors, and if 
we do not get the short distance behavior right, 
it is not likely we will get the interpolation right. 

Thus we either must abandon the strong cou- 
pling expansion or be optimists and look for a 
version for which summation and analytic con- 
tinuation produces the correct result as g ^ 0. If 
we can represent it as a string theory, an appro- 
priate treatment of that theory may well do the 
analytic continuation for us. A solid result which 
speaks against this possibility is the roughening 
transition - I refer to for a discussion of this 
but note only that the real cause of this was the 
lattice discretization, and that a continuum string 
theory will not have it. 

The most direct and explicit version of this idea 
would be to develop a continuum strong coupling 
expansion. In Z) > 2 we need a regulator, but in 
principle any continuum regulator which allows 
us to make sense of / DA without the Boltzmann 
weight might be used. Although it would depend 
on the choice of regulator, it would surely be a 
more natural expansion than the lattice expan- 
sion, and conceivably free of spurious phase tran- 
sitions. This is one way to make the concept of 
'QCD string' precise. We will make contact with 



3 



my original definition if we can associate terms 
with surfaces, if the weight for each surface is de- 
termined by local rules, and if we can reproduce 
the weight with a local world-sheet action. 

Such a string theory has been developed by D. 
Gross and W. Taylor for two-dimensional Yang- 
Mills theory. |Q They gave local rules for the 
weights and it can also be defined by a continuum 
world-sheet action, as described by G. Moore in 
his lectures. | |Tc| , |2^ 

There is no obvious reason that a similar ex- 
pansion for a regulated higher dimensional theory 
could not exist. The regulator would have to act 
solely on the measure and preserve gauge invari- 
ance - thus the likely candidates are stochastic 
regularizations. I feel this possibility is the 

most important lesson to be drawn from the work. 

Will it be QCD? In other words, will the g ^ 
limit reproduce weak coupling physics? It is dif- 
ficult to get at this question using the expansion 
itself - a priori, even if it converges, we will need 
more and more terms to get a result as g ^ 0. 
We expect the world-sheet theory to be strongly 
coupled at short distances, jl^ Conceivably, we 
might end up with an exact string reformulation 
of QCD, in which (say) the one-gluon exchange 
graph is uncalculable. 

A better approach is to use other methods 
to study the analyticity of the partition func- 
tion. Combining this information with the known 
finite radius of convergence of the strong cou- 
pling expansion will show its validity. Clearly 
we should start by answering this question in the 
two-dimensional case, which we can solve exactly. 

2. Tw^o dimensional Yang-Mills theory 

There are many ways to solve the two dimen- 
sional theory, and it is worth doing it in more 
than one way, first because the techniques come 
up elsewhere, and second because it is conceivable 
that one of them will provide insight for higher 
dimensions. 

YM2 is very easy to work with, first because 
the action / (Px tr only depends on the vol- 
ume form, and second, as pointed out long ago 
by A. Migdal, there is a lattice definition which is 
equivalent to the continuum theory. The Boltz- 



mann weight ('heat kernel action') Z{U,A) for a 
plaquette of area A and boundary holonomy U 
is just the continuum YM2 path integral on the 
plaquette. Let the plaquette be the disk with co- 
ordinates {r,9) and 7rr^ < A; the A dependence 
can be determined as Hamiltonian evolution with 
r as time. Fixing Ar = gauge and imposing 
Gauss' law, the wave function depends only on 
the holonomy U ~ expi J dO^'^Ag. The Hamilto- 
nian is 




This is the only appearance of and from now 
on we normally choose units of length in which 
g^ = 1. E" ^ tr f'Ud/dU generates left rota- 
tions of U and represents the Lie algebra u{N). 
Thus acting on a wave function which could be 
any matrix element of an irreducible representa- 
tion R, il){U) = Df\u), H = C2(i?), the second 
Casimir (normalized so that C2 (□) = N) . Gauge 
invariant wave functions ipiU) — ipigUg^^) are 
class functions, which can be expanded in char- 
acters XBiU) = tr D(^HU). Let xniU) be 
the character of the irreducible representation R 
with X-r(I) = dimi?; the measure J DU is nor- 
malized so / DU XRiU)xsiU+) = Srs- These 
results combine to determine the path integral 
Z{Ui, U2] A) on a cyhnder of area A and bound- 
ary holonomies Ui and U2'- the final ingredi- 
ent is the boundary condition Z{Ui,U2;0) = 
SiUi,U2) = J2rXr{Ui)xr{U:^) on class func- 
tions, so 

Z{Ui,U2;A) - ^ Xfl(t/i)xfl(C/2+)e"*^''^^''^-(5) 
Red 

{G is the set of irreducible representations of 
G). We can then set U2 — I, the holonomy 
for a zero area plaquette, to get Z{U;A). By 
writing xr(UV) = tr D^^'> {U)D^'^'> (V) and us- 
ing / dC/^(f(C/)7^(f([/+) = ^SnsS^Aj one 
can check the 'self-reproducing' property 

J dVZ{UV-^;Ai)Z{VW;A2) 

= Z{UW;Ai+A2). (6) 
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This implies that the lattice integral is invariant 
under subdividing the discretization and thus this 
is equivalent to the continuum limit. 

Perhaps the simplest result is the calculation of 
the partition function on a Riemann surface due 
to A. Migdal and B. Rusakov. [Q We can form 
a genus g surface by identifying the edges of a 
4g-gon as AiBiA^^B:[^A2B2A^^B2^ ... and it 
is an entertaining calculation to check that 

Zg(A) = ^ (dimi?)2-2g e-^'^^^^'. (7) 
Red 

We specialize this to the group U{N) to study 
the large N limit. This choice turns out a bit sim- 
pler than SU (N) but one might worry that we are 
introducing extra degrees of freedom which will 
confuse matters. The C/(l) factor is a gauge the- 
ory with gauge coupling /N —> in the limit, so 
perturbatively this completely decouples. How- 
ever, U{N) is not a direct product but rather 
a quotient by Z^v: one has the exact sequence 
Zn ^ SU{N) X U{1) ^ U{N) where the first 
map is n ^g2Trm/JV-|^^ g-2Trm/JV^ sccond is 

(e*^, U) e^^U. In (|^) this means the sum over 
representations is projected to those with U{1) 
charge equal to the A^-ality of the SU{N) rep- 
resentation, modulo N. This does affect exact 
results but in a qualitatively unimportant way. 
(It turns out that the g = result is unaffected; 
this will follow from the saddle point calculation 
below and can also be understood from the string 
point of view [sil]). 

To get a feeling for it, let us look at the first 
terms of the series: 

Zg = l + 2iV2-2gg-4 

+2(i7V(iV + l))2-2ge-4(2+2/JV) (8) 
-f2(iiV(A^- l))2-2ge-f(2-2/w) 

+ (iV2 _ i)2-2gg-42 ^ ^ ^ ^ 

Restoring g^, this is an expansion in e^p—g^A, 
a parameter which is zero in the strong coupling 
limit. Thus this is a strong coupling expansion, 
and we can study its validity as g 0. This 
is the limit ^ ^ 1, and this extrapolation 
already sounds more attainable than 1/g^ oo. 
This type of improvement is common in character 
expansions. 



As was explained here by G. Moore, the work of 
Gross and Taylor |24) shows us how to associate 
the 7V2~2Gg-nA/2 ^g]-j32g jj^ f^Y^Q expansion of 

(P) with n-fold branched covers of the original 
surface by genus G world-sheets. Let us look at 
g = and G = 0: the 0{e~^/^) terms are easily 
associated with single covers, while the 0{e~^) 
terms in Z are 

{2N'^ + {^A^ -2A-l)N^ + 0{N^))e-'^ (9) 

The first term is the disconnected piece from the 
e-^/2 diagrams. The 0{N^) terms must be re- 
produced by double covers of the sphere, and the 
rather non-trivial weight can only be reproduced 
by defining two types of branch points, which we 
can call 'ordinary' (coming with a power of A) 
and Ti-points' (with weight 1). This suffices at 
higher orders, and one consequence of this is that 
the 0(e~""^/2) term associated with n-fold covers 
is weighted with a polynomial in A with order the 
number of branch points, 2n — 2. 

The expansion has a very different character 
for g = 0, g = 1 and g > 1. At fixed order in 
1/iV, for g > 1 we have finitely many terms, and 
no subtleties arise. For g = 1 one can express the 
results in terms of modular forms (ij,^^ and the 
expansion is valid a,s g^A —> 0. 

The situation at g = is not immediately clear. 
There will be an O(N^) free energy, as is true in 
higher dimensions. This is the term which we will 
relate to genus zero string world-sheets in a string 
interpretation. We thus expect it to be 'classical' 
in the same sense that tree diagrams in quantum 
field theory are calculable using classical field the- 
ory. This is a central observation and underlies 
many of the other formalisms which have been 
proposed to solve large N field theory; we will 
come back to it again. 

Let us now reformulate (Q) in a way which will 
provide exact results. We just need to make the 
group representation theory explicit, and rather 
than do this abstractly, let us find the eigenfunc- 
tions of tr directly, i.e. solve the quantum 
mechanics defined by (jj), following (This 
is also a standard mathematical approach, as in 

up.) 

Since we are most interested in class functions. 
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we change coordinates on the group manifold to 
Uij = gikZkg^j ■ The invariant volume element in 
these variables is 

\/X= |A(z)|2 = A(z)2 (10) 
where /^{z) = Wi^jizi - Zj) and K{z) ^ 

n,<,sin^ = A(z)/n.^f"''^'; The "ra- 
dial" components of the metric are simply hij = 
6ij . Thus on wave functions independent of g 



^ A2 dB^ dOi 

i 

We can rewrite this as 



H 



E 



1— A-- 



d^A 



A 



(11) 



(12) 



The second term, after some calculation, is found 
to equal -N{N'^ - 1)/12. 

Thus, we can redefine the wave functions by 



i' = Ax, 



(13) 



and arrive at a theory of N free fermions on the 
circle. The boundary conditions are also deter- 
mined by this redefinition; they become periodic 
(antiperiodic, respectively) for N odd (even). An 
orthonormal basis for wave functions is Slater de- 
terminants 



= det z"^ 



(14) 



with energy E = Y.i^l " ^(^^ - 1)/12- The 
ground state has fermions distributed symmetri- 
cally about n = 0, and energy zero, so the Fermi 
level np^iN ~ l)/2. 

Going back to the original wave functions, we 
have rederived the Weyl character formula: 



Xn{z) 



deti<tj<jv z"^ 

1 , 7 — 1 — n F ' 

deti<ij<Ar zl 



(15) 



In terms of roots and weights, the indices rii 
with ni > > . . . > njv are the compo- 
nents of the highest weight vector shifted by half 
the sum of the positive roots (usually denoted 
IJL + p) where the basis of the Cartan subalgebra 
is just {Hi)jk = SijSjk- In the language of Young 
tableaux, if hi is the number of boxes in the z'th 
row, = (TV - l)/2 + 1 - i + ft-i. 



The U{1) charge is Q — J2i We can change 
this by a multiple of N by shifting all the fermions 
ni ^ rii + a, but Q mod A'^ is correlated with the 
conjugacy class of the SU {N) representation. 

We will do more with this later, but for now 
we simply adopt the labelling scheme {n} for 
representations, the formula C2 = '^n^, and fi- 
nally the dimension formula for a representation, 
computed by rewriting the determinants as Van- 
dermondes and using rHopital's rule to take the 
limit: 



dimi? = ^hnx«(C/) =n 



i>j 



^- 3 



(16) 



Substituting into (^ we find 



N 



lim Zg(A) = cxp[iV2i^(A)] = 



En 

m^iij i>j 



A 



exp - 



N 

2N ^ 

i=l 



(17) 



As pointed out by Rusakov [Q , for g = this 
expression is such that the sum over {rii} can be 
determined by saddle point: it is 



e^p[N^F{A)] = lim 



N- 



exp- 



-N^S^ffin) (18) 



with 



Seffin) = -^E^og 



A^ 



(19) 



A 
2N 



Thus our general expectations that 'leading large 
A"' = 'genus zero string theory' — 'classical the- 
ory' are confirmed. 

Seff contains a repulsive or 'entropic' term 
which favors non-trivial representations and the 
sum is exactly the same as a hermitian one-matrix 
integral with the eigenvalues replaced by quan- 
tized variables rii/N . Since the quantization is in 
units of 1/A^ one might expect it not to affect the 
leading order saddle point calculation of F, and 
this is almost true. We thus introduce a scaled 
variable p — n/N and spectral density 

1 



(20) 
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Figure 2. Spectral density on S*^ in both phases. 



in terms of which 



Seff{p) 



dpdq p[p)p{q) log \p - q\ 
A f 



(21) 



This is the effective action for the simplest of 
matrix models, the gaussian integral, which we 
could just do directly: 



/ 27r 

Vna 



(22) 



This is a very simple answer, but it appears to 
have no relation to the original sum (0)! The 
original sum over representations or the string 
representation derived from it seems to have been 
very misleading. Furthermore, there are indepen- 
dent arguments for the simple answer (22), such 
as the suppression of non-trivial contributions in 
the heat kernel on the group manifold evaluated 
at time t — A/2N or of instanton contribu- 
tions in a direct evaluation of the path integral 

The resolution of this paradox is that we ne- 
glected a crucial consequence of the discreteness 
of the variables Ui |18|: the bound 

p{p) < 1. (23) 

The maximum density is attained if the rii are 
successive integers. 

The saddle point for (E|) satisfies 



A 

2P = 



dq 



■piq) 



(24) 



and, 



Ignoring 
A 
2^ 

For A > A^ 



is the semicircle 



pip) = 7rv/4M-p'- 



(25) 

^ this violates the bound and 



we must find a new saddle point, enforcing the 
constraint by hand. Thus, letting p{p) = 1 for 
\p\ < b and integrating this range of p in ( |24|) , 



A . p-b 

-^P - log — r 

2 p + b 



dq 
p-q- 



p{q). (26) 



The general solution of such linear inhomoge- 
nous equations is known and the result is given in 
|l8|] . It is expressed in terms of elliptic functions, 
and its series expansion reproduces (||). A graph 
of pij)) in the two phases is in figure 2. 

For g^A < tt^, the strong coupling expansion 
fails. The partition function is non- analytic in a 
particularly drastic way - its two branches are 
completely unrelated, because 'turning on' the 
constraint is a non-analytic operation. This 'large 
N transition' is a consequence of taking the limit 
N oo and does not have a direct analog at 
finite N, as will be clearer below. 

The original result of this type was due to Gross 
and Witten and was generalized by Brezin 
and Gross H] to the matrix integral 



Z{A) 



dUe 



Nf3ti- AU+A+U^ 



(27) 



This integral is the generating function for the 
link integrals such as (||) used in deriving the orig- 
inal strong coupling expansion and the conclusion 
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was drawn that this expansion would be invalid at 
small g. Now we have learned that the transition 
is not a lattice artifact. 

The double scaled theory around the new tran- 
sition is the same as that for the Gross- Witten 
transition, with the roles of weak and strong cou- 
pling interchanged. [p3| This may be intuitively 
plausible from figure 2. 

In terms of the strong coupling expansion, the 
YM2 transition is signalled by non-analyticity 
from summing the polynomial prefactors in (||), 
and thus there is a string theory explanation of 
the transition |51|: the sum over the number of 
'ordinary' branch points diverges at the transi- 
tion. 

To make a similar argument in higher dimen- 
sions, we would need to be able to control the 
signs which appear in the expansion. This is 
probably not realistic as it is known that the signs 
must drastically change the asymptotic behavior 
of the string partition function to be subexponen- 
tial in the area. |^Jlll 

All this is a rather serious blow to the QCD 
string idea as there is no other convincing argu- 
ment for an exact string reformulation. Thus it 
is essential to understand whether this affects the 
physics, whether it persists in D > 2, and if so 
what we can do about it. To say anything about 
these questions, we need to translate our results 
to a formalism which could work in principle in 
D>2. 

3. Collective Fields and Bosonization 

The idea that 'the large N limit is classical' can 
be formalized in many ways. At the heart of it is 
factorization: 

(tr Oi tr O2) = (tr Oi)(tr O2) + 0{N°). (28) 

This follows directly from the topological expan- 
sion - the factorized term is a disconnected di- 
agram of 0{N'^). It is as if the functional inte- 
gral was dominated by a single saddle point - but 
one should realize that this is not literally true 
in £) > 2, and that only C/(A^)-invariant quanti- 
ties have well-defined expectation values. Neither 
has it been proven beyond perturbation theory in 
D > 2, but its truth in both weak and strong cou- 



pling expansions is a rather convincing argument. 

We thus would like to formulate the large 
N limit as a classical theory whose configura- 
tion space is parameterized by the expectations 
tr Oi and whose ground state is the solution of 
some 'equation of motion.' One candidate for 
this equation is the factorized Schwinger-Dyson 
or Migdal-Makeenko equation, as A. Migdal ex- 
plained in his lectures: in the continuum. 



^ J dy'' 5{x ~ y) WiL^y)W{Ly^) 



(29) 



This equation is quite tractable in two dimen- 
sions because of the area-preserving diffeomor- 
phism invariance. For present purposes, however, 
it is easier to work in a canonical formulation, 
with which we can make contact with the results 
of the previous section. 

There is a general procedure for finding a 'clas- 
sical' Hamiltonian and phase space reproduc- 
ing the large N limit of a general field theory, 
the collective field theory of A. Jevicki and B. 
Sakita. The most general and complete treat- 
ment of the canonical formalism is due to L. Yaffe 
and collaborators, and I highly recommend [p6|,p[ 
for the reader's further study. (We unfortunately 
did not have time in the lectures to do justice 
to this.) The formalism applies in particular 
to gauge theory in any dimension, and although 
there are still not many analytic results from it 
in D > 2, there is a good deal of numerical evi- 
dence that it makes sense and properly describes 
the large N limit of the regulated theory, (e.g. 
||^,|3^ and references there.) 

Before plunging into details, let us repeat our 
primary question: the strong coupling expansion 
provides a fairly explicit description for a candi- 
date ground state of our gauge theory, so why not 
just use it? 

The simplest derivation of the collective Hamil- 
tonian is to change variables in the quantum 
Hamiltonian to invariants. In Z? = 1 -|~ 1 this is 
(|T]]), H = — ^diVhh''^ dj and the invariants are 
Wn — tr [/", or the associated spectral density 
Wn^ Jd9 p(6')e™^ 

If we want to write the Hamiltonian using 
canonically conjugate variables, we need to make 
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a wave function redefinition like (|f3|). This 
can be seen by the foUowing heuristic argument: 
we want the invariant inner product (xlV') = 
/ dpVh/x*{p)i^{p) to be simphfied by the redef- 
inition h'^/'^ip ip, so that the self-adjoint mo- 
mentum operator 11 = h'~^/'^{5 / 5p)h'^^^ — > 5/5p 
which has canonical commutation relations with 
p. Carrying this out produces an 'effective poten- 
tial,' which can be expressed in terms of invari- 
ants. 

nj4£ 



Let us quote the by now standard result 1 34 



He 



1 

2N 



d9 pidelV) 



(30) 



where {Tl{e), piO')} = 6(6-9'). p{9) is a spectral 
density and thus there is a constraint J p — N 
(which could be implemented with a Lagrange 
multiplier) as well as the inequality p > 0. For 
now, we will not redefine p — > p/N, so we can 
better explain the relation to the topological ex- 
pansion as well as to the original quantum theory. 

A simpler set of variables are the chi- 
ral combinations p± — dgll ± T:p satisfying 
{p±{e),p±i0')} = ±27TdeSi9-0'): 



H = ( de p\- 

UttN J ^+ 



(31) 



Classical time evolution under this Hamilto- 
nian reproduces the large N limit of real time 
quantum evolution in group quantum mechan- 
ics. The equations of motion derived from ( |30| ) 
are Euler's equations for a one-dimensional fluid. 
(The velocity is v = dglT). In the variables p± 
they decouple: 



dtP± + -^P±dep± = 0. 



(32) 



This is often referred to as the Hopf equation. 
Our YM2 problems are defined in two 'Euclidean' 
dimensions, and we will have to take this dif- 
ference into account below: we will see that we 
should take t — it and 11 — > —ill. 

In £> = 1 + 1 one can derive the collective 
field theory from the simpler free fermion solu- 
tion. (Complete details are given in so we 
only give a summary in these notes.) The invari- 
ants Wn are expressed in terms of the fermions 



as tr [/" = "^iZ". We second quantize, in- 
troducing the creation and annihilation opera- 
tors — Sm.n and non-relativistic fermi 
fields ipiz) = Y.n ^nz"-. We then have tr [/" = 
/ dzz-"-'^^p+{z)4>{z). Thus p(z) = ^+(z)-!/;(z) 
is the spectral density and we would like to talk 
about it as a classical field. 

The general answer to this type of problem 
is bosonization. However, the classic Coleman- 
Luther-Mandelstam bosonization can only be ap- 
plied to a relativistic fermion. We can argue as 
follows in the large A'^ limit: if we agree not to 
consider operators tr with |n| ~ we can 
regard excitations of the two fermi surfaces as 
completely independent, and extend the fermi sea 
below np = {N — l)/2 to n = —00 and similarly 
extend the sea above —np to n = 00. The result- 
ing system is kinematically a complex relativistic 
fermion. (This is a standard argument in many- 
body theory.) The ability to decompose the ex- 
citations into right movers ip{z) and left movers 
'(/'(z), describing excitations from the two fermi 
surfaces, corresponds in the original group the- 
ory to decomposing U{N) representations into a 
tensor product of a 'chiral' representation built 
from finite tensor product of fundamentals with 
an 'antichiral' representation built from the anti- 
fundamental. Rewriting the Hamiltonian 



Hnr = J dz z'^d^j^dip 

in terms of the relativistic fields produces 



(33) 



H = §dzl 
+ §dzl: 



ip+dx/: : +j^z^ : dip+d^j : (34) 



where Lq and Lq are the standard conformal 
field theory Hamiltonians. Applying relativistic 
bosonization —: e*"^ : produces 

H^j i p</,)2 + (9^)2] + _L [(90)3 _(a^)3]. (35) 

This is (HI) with p+ = N/2 + dcf) and p_ = 
—N/2 + dcj). (We are taking as a standard free 
field and thus implicitly defining the time deriva- 
tives in d = ^{do — di) using free time evolution 
generated by io + -^o-) We see that it is valid to 



all orders in 1 /N (if we attend to quantum normal 
ordering) . 

Let us explain the connection between this for- 
malism and gauge strings. The strings are the 
'quantum' fluctuations of the field </>, of 0(1). 
Since tr [/" = q;„ + a_„ after bosonization, we 
see that we reproduce the picture from the strong 
coupling expansion if we identify the excitations 
of the n'th left-moving bosonic mode with n- 
winding strings, and the n'th right-moving mode 
with — n-winding strings. The leading 0{N^) 
part of the Hamiltonian ( ^ ) preserves string 
number and gives a string energy proportional to 
the absolute value of its winding number, while 
the 0{1/N) piece is a three-string interaction. 

For present purposes, we can identify world- 
sheet and target space time, and think of the 
world-sheet Hamiltonian as / day^X'{ay . The 
complete theory is an interacting 'string field the- 
ory.' 1 41 1 Factors coming from [a„, = n make 
the interaction amplitude for an n winding string 
proportional to n. This can be reproduced by 
a simple vertex which splits or joins strings with 
equal amplitude at each target space point. These 
are the 'ordinary' branch points in the covariant 
language. 

What makes the discussion intuitively straight- 
forward is the identification of world-sheet and 
target space coordinates, making this a tempt- 
ing assumption in Z? > 2 as well. However, the 
canonical approach is not the best way to derive 
a string because it loses target space symmetry. 
A major advantage of a string reformulation is 
that what would have been a sum of diagrams in- 
volving any number of vertices, becomes a single 
diagram in a covariant approach. 

The canonical formulation does allow us to eas- 
ily compare the string approach with exact re- 
sults. Let us focus on the question: does the 
genus zero string theory correctly reproduce all 
'classical' (leading large N) results? These are 
0{N) terms in Wilson loop expectations and thus 
in (/) - compared to this, the commutators are sub- 
leading and thus 'quantum.' 

If we perturb a coupling by 0{N) (or add 
such a source), the resulting perturbation of the 
ground state (j) is also 0{N) - we have changed 
the 'classical' theory. Now there is no a priori 




Figure 3. A lattice with one plaquette at each 
time. 



reason why such results cannot be reproduced by 
string theory, even though the assumption that 
the amplitude of the perturbation is 0(1) is now 
false. However there may be 'non-perturbative' 
structure in the theory, in other words structure 
which appears only for 0{N) shifts of the fields, 
and we will miss it in the string theory. 

A model which can be solved either by min- 
imizing H or by 'string perturbation theory' is 
the following, studied by Wadia: pj] 



H 



N 



-tr 



N 



(36) 



The potential term is a simple approximation 
to the space-like ('magnetic') plaquettes oi D > 2 
lattice gauge theory - we could get this model by 
starting with the lattice in figure 3, and taking the 
lattice spacing in the time direction to zero. The 
feature of D > 2 gauge theory we hope to capture 
with this is the following: at very weak coupling, 
the gauge field U{x) — expiA{x) is roughly Gaus- 
sian, and each mode A{k) has a Gaussian wave 
function with width g^jfeja. The potential pro- 
duces a similar behavior here and we can think of 
this as the dynamics of a single mode with \k\ ^ a. 

To find the ground state, we can sim- 
ply minimize ( pO[ ) with a potential term 
/ d9 p{9) (^1^(1 — cobO) — fi) to find the ground 
state 



N 



1 



(1 



(37) 
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with determined by J p — N. Again we need to 
enforce the constraint p{6) > on this solution by 
hand and there are two cases: for (?^a^ > it 
is never saturated, while for g^a^ < 4/0? there 
is a region around 9 = tt where it is. Thus this 
model has a large A'' transition as well. 

We could also have expanded H/g"^ in /? = 
l/g^a"^ to get the Hamiltonian strong coupling 
expansion, and resummed this expansion. In the 
string picture this is a source and one picture we 
can make is that we are summing world-sheets 
with k 'holes' at order l/{ga)'^'^, and integrating 
over the time in the target for each hole. Since 
the source is 0{N), terms involving n cubic string 
interactions and n + I sources are leading order. 

To preserve the D > 2 analogy, we can also 
postulate continuous world-sheet embedded in 
the 'magnetic plaquette.' It will have different 
weights from those in the original expansion: es- 
sentially, there are no branch points on the new 
world-sheet. One could enhance the analogy even 
more by taking the heat kernel action for the 
magnetic plaquettes, so the world-sheets would 
be generated with the same weights. This does 
not change the qualitative behavior. 

Summing all disconnected diagrams will pro- 
duce a coherent state for the ground state, of the 
form 

= exp ^a_„ + -^a_„ j |0) (38) 

with Wn determined from p{9) in (p^). 

The exact result is analytic in /3 near 0, so sum- 
ming the string diagrams will reproduce it. If we 
try to go past the transition by analytic contin- 
uation, we will produce a complex p, the con- 
tinuation of defined by simply ignoring the 
constraints. 

Was this a failure of string theory or of the 
strong coupling expansion? Really, it is a fail- 
ure of both. Built into the free string theory is 
the assumption that we can vary the occupation 
number of each winding number of string inde- 
pendently. What we see is that this can fail if the 
occupation number or amplitude is 0{N). This 
is the 'non-perturbative structure' we alluded to 
above. We can trace the failure back to the step 



where we decoupled the two fermi surfaces, a nec- 
essary step in deriving the bosonic theory. 

Now the equivalence of ( |30| ) and (|33| ) pro- 
vides a more general rewriting the theory in 
terms of bosonic variables. We could call it 
'non-relativistic bosonization.' This is particu- 
larly interesting in the present case, since the ar- 
gument involving decoupling the fermi surfaces 
broke down at the turning point (where p just 
reaches zero). This bosonization is exact in the 
classical limit, and (with careful treatment of the 
turning point) quantum corrections can be de- 
rived from ( pO| ) as well. This is the case relevant 
for the c = 1 matrix model. p3| | 

We still cannot use this to fix up our original 
string theory, because we still need the constraint 
p{9) > which cannot be expressed in string lan- 
guage. What it does suggest is that a broader idea 
of string theory might exist, and we will make 
some comments about this below. 

In all the problems we discussed so far, and 
the one-matrix integrals, large N transitions all 
appear exactly where we begin to saturate a con- 
straint, p{9) > or p{p) < 1. Let us see why these 
are two forms of the same underlying constraint. 

Since the eigenvalues of U are classical fermions 
in the large N limit, we should be able to spec- 
ify their positions and momenta simultaneously. 
Since they are indistinguishable, the state of the 
total system is given by a phase space density 

p(0,p) -^<5(0-0,)'5(p-pO- (39) 

i 

Given this, p{9) = J dp p{9,p) and p{p) of (^ ) 
is jj J d9 p{9,p). The fermionic nature of the 
eigenvalues is exactly captured by 

N 

< p{9,p) < —. (40) 

Integrating d9, we see that the compactness of 9 
is responsible for the upper bound on p{p). 

The phase space description is more general 
and in many ways simpler than the collective field 
theory. | |4^ It can be derived from the quantum 
(non-relativistic) fermi theory using 

p{9,p)^ [ dae'^P ^+{9-a/2N)^{9+a/2N).{41) 
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Any one body operator can be written directly as 

f dedp f{e,,p,)p{e,p). (42) 
i •' 

(This can be used quantum mechanically and pro- 
duces Weyl ordered operators.) 

For the operators, the large N limit is the stan- 
dard classical limit. This is not necessarily true 
for the states - assembling N classical fermions 
will produce states in which p at each point is ei- 
ther or pmax = N/2iT^ but in general one can 
construct states in which < p < Pmax- An 
example is the state corresponding to a specific 
representation R, i.e. with wave function det ' . 
Such states are not coherent states of the form 
( p8| ) and do not usually come up in practice. 

The commutator of operators becomes Pois- 
son bracket, a free particle bracket {/, .g} = 
2Tr{dpfdgg — dgfdpg). Time evolution is p = 
{H,p}. 

Collective field theory is derived by assuming a 
form 

p{9,p) = ll iorp^{0)<p<p+{e) (43) 
1^ otherwise. 

For example, H — J dp p{9,p)\p'^ reproduces 
(^l|). Again, the simple form of the state is an 
assumption, which is easily seen to be correct for 
problems such as ( |36| ) . Several generalizations are 
possible. One can have a fermi surface of arbi- 
trary shape, and a collective field theory using 
several functions Pi{0) to describe it. There is 
also nothing special about the parameterization 
p{9) and one could use 9{p) or [9,p) depending 
on a parameter s to describe the surface. 

Free fermions are special to D — 1 + 1 (and this 
system), but a language which in principle could 
generalize some of this to I? > 1 + 1 is known. 
We can construct the phase space as the orbit of 
a reference state under the action of a 'coherence 
group,' obtained by exponentiating the Lie alge- 
bra of observables. Thus the entire kinematics of 
the limiting large TV theory is contained in the 
geometry of this group. 

The discussion of started from the ba- 
sic observables of collective field theory, p{9) de- 
rived from the invariant operators tr [/" and 
pdXl{9) from ifr [EV' + U'^E). These generated 



a semi-direct product of a Virasoro algebra with 
the modes of p{9) and the corresponding group 
Diff(S'^) semi-direct product the additive group 
of functions p{9). The algebra generalizes directly 
to £) > 1 + 1 and in the next section we will be 
able to say something about the coherence group. 

More recent work on these lines uses 
all the one-body operators (^) to generate the 
group of symplectic diffeomorphisms on phase 
space, "SDiff(2)." The idea that the canonical 
formalism is simpler if we consider all invariant 
operators made from U and i?, not just low pow- 
ers of E, has not been studied in D > 2, and 
might be valuable there as well. 

So far, we have seen a relation in _D = 1 and 
D — 1 -\- 1 between large TV transitions and con- 
straints on the invariant observables. Similar con- 
straints exist in higher dimensions. Essentially 
nothing is known about the phase space version, 
but the analog of p{9) > is clear. Given any loop 
L in space-time (or lattice), we can consider its 
holonomy U{L) and associated spectral density 
Pl{0)- this should satisfy the same constraint. 
We will look at this in more detail in the next 
section. 

It is natural to conjecture that large N transi- 
tions are always associated with a change in the 
application of the constraints. Now, to have any 
hope of making a definite statement in higher di- 
mensions, we need to find some qualitative con- 
ditions for the transition, that do not require ex- 
act results. Thus we might hypothesize that if 
for coupling gi the ground state saturates a con- 
straint, and for 172 it does not, there is necessarily 
a large N transition in between. 

This is an attractive hypothesis because it can 
be checked at strong and weak coupling where 
we have some control over the higher dimensional 
gauge theory. For strong bare coupling, pLiO) 
(for a loop on any scale) should saturate no con- 
straints {p > V0), while at weak coupling and 
short distances the gauge field fluctuations are 
highly suppressed and we expect pL for small L 
to have support sharply peaked about the origin 
(with width ~ g'^a/\L\), no matter what regula- 
tor we use. This has been seen in the numerical 
studies |3^J^ and the finite radius of convergence 
of both diagram expansions makes it accessible 
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n-m 



m 



analytically. The only doubt we might have re- 
gards the weak coupling phase, where to be pre- 
cise we must be able to distinguish p{9) = from 
a very small but nonzero value, but let us imagine 
we can prove that (say) p{-k) ~ for small loops. 

One might think (and at Trieste I advocated 
the idea) that this would prove the existence of 
the large A'^ transition in higher dimensions, and 
the failure of the strong coupling expansion. We 
might expect for small loops to become com- 
plex, like (37), which looks very unphysical. 

However, there are subtleties in this, as one 
realizes on considering the Wilson loop expecta- 
tions in two dimensions on the plane. These are 
determined directly from the Boltzmann weight, 
Z{U] A) = Z{U, 1; A) as in §): 

W„ = I dU tr C/" ^Xfl(C^)dimi?e-^^2(^) 

•' R 

n 

= ^(-l)™dimii;e-^^^(^) 

m=0 



The sum is over the representations in the 
above figure. (This can be derived from fermion- 
ization: tr C/"|0) = &V2-n+™^-i/2-™|0).) 
Since the sum is finite, this will be analytic in 
A. (The string diagrams can have up to n — 1 
'ordinary' branch points). If we compute p{9\ A) 
from these, we will find that it satisfies the con- 
ditions of the hypothesis: for g^A < 4, p{6) has 
a gap around 9 — tt, while for g^A > 4 it does 
not. Il^] We might call this a 'pseudo-transition'. 
It is visible in observables, but we need to look 
at Wilson loops winding arbitrarily many times 
to see it: at fixed g'^A < 4, the W„ 



,-3/2 



We thus see that non-analyticity in the cou- 
pling for individual Wilson loop expectations and 
the existence of phase transitions cannot be con- 
cluded solely from the behavior of the spectral 
density, as pointed out long ago by B. Durhuus 
and P. Olesen. Actually this does not directly 
contradict the hypothesis as stated, because the 
states in question are not ground states of the 
Hamiltonian. Nevertheless our belief in it may 
be somewhat shaken. 

One can compute p{9; A) more directly from 
the collective field theory. The problem was ac- 
tually solved this way in [|l9|, without ever men- 
tioning collective field theory - they derived ( ^ ) 
directly from the Migdal-Makeenko equation. 

Let us do 'radial quantization' on the plane, 
embedding the Wilson loops at fixed radius r, 
and and then go to the canonical formalism us- 
ing A = 7rr^ as time. Let us also finally redefine 
p p/N from (|0|). p{9; A) is thus the result of 
evolution for Euclidean time A from the bound- 
ary condition p{9;0) = 5{9). It is simplest to 
think of going to Euclidean time in the action 
S = J dtUp — H, by taking t = iA. To preserve 
the commutation relations in the Hamiltonian we 
should restore the first term by taking H —ill. 
It is then convenient to redefine 



— dgll ± inp 



(44) 



The equations of motion are the same (except for 
the 1/iV): 



dAP± +P±dgp± = 0. 



(45) 



as 



oo, rather than falling exponentially. 



A nice simplification is now possible: since H 
and p are real, the two equations are complex 
conjugate, and and if p+(9) is analytic, we have 
reduced the equation to a first-order ODE. Since 
p{9) and H(0) are only prescribed on the real axis, 
there is generally no difficulty in finding such an 
analytic function. 

We don't know H at A = - it is not deter- 
mined by the quantum initial condition i){U) = 
5{U,1). In general, the classical limit of a Eu- 
clidean time quantum problem becomes a bound- 
ary value problem, where we specify (say) initial 
and final positions. The initial and final momenta 
are then determined only by solving the equa- 
tions of motion. In general, more than one solu- 
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tion can exist with the same boundary conditions, 
and this is a possible source of phase transitions. 
All this applies to the large N limit and in the 
present case we must specify the limiting behav- 
ior of p as A oo: it is the strong coupling 
vacuum p{9) — l/27r. (Similar considerations ap- 
ply to the large N limit of the Itzykson-Zuber 
integral, which as shown by A. Matytsin, can be 
expressed in terms of the same boundary value 
problem. 

Given only the boundary conditions, the ex- 
plicit solution is not easy to find. However there 
is one more fact about the solution, which greatly 
simplifies the problem. After extending p± to be 
analytic functions in the complex 9 plane, they 
fall off as p± ^ 1/9 for large 9, much like the re- 
solvent tr (z — C/)^^. Why this should be is not 
at all obvious from what we have said so far, but 
this assumption makes the problem easy, because 
then n at A = is determined by a dispersion re- 
lation, and we have an initial value problem. We 
can check the assumption at the end by verifying 
that it produces the correct A — > oo limit. 

The final ingredient is to take periodicity in 9 
into account, which we can do by periodizing the 
initial condition. The analytic function satisfying 
Im/ = nJ2n ^{0 ^ 2nn) is then 

/W = ^cot^. (46) 

The general solution of the Hopf equation 
can be derived by first representing the bound- 
ary p+{9) in phase space parametrically as 
(p+(^), 9{£,)), implementing free motion 6'(^; A) — 
9{^; 0) -I- Ap+{^), and solving for ^. Although the 
pictures are simplest for real p+, the solution is 
perfectly valid for complex p+ . 

Thus, given an initial condition p_|_(0; A = 0) = 
f{9), the solution is 

p+{9;A)^fmA)) (47) 
where ^ is the solution of the implicit equation 



(48) 



We can then recover p = ^Imp_|_ — — 

This cannot be solved in closed form, but the 
comparison of p{9) with (jij) for A ^ oo, and the 
verification of the properties we described are not 



hard (and done explicitly in |19|] ). For example, 
to check that p{9; A) 1/2tt as A oo, we can 
take the imaginary part of (^) and check that 
Im^ -A/Att. 

What is the essential difference between this 
problem and the sphere partition function? In 
this language, the sphere problem is formally very 
similar: it is a two-point boundary problem where 
both initial and final boundary conditions are 
p{9) = N6{9). 

Again, boundary value problems are not easy 
to solve, so we ask if another trick can convert 
this into an initial value problem. We can al- 
ways extend p± into the complex 0-plane as an- 
alytic functions, but now they cannot fall off as 
P± simply because 11 is different in the 

two problems. I have not found such a trick 
and since I just want to illustrate the qualita- 
tive aspects of the problem, let us instead solve 
the problem for fermions on the real line. In 
other words we replace YM2 by hermitian ma- 
trix quantum mechanics. In this case we already 
know the solution of the quantum problem with 
boundary conditions tp{M;t = 0) = S{M): it is 
\t) oc exp — iVtr AP/t. Using this we can compute 



Pt(A) = (T-i|i^5]<5(A-AO|i> 



(49) 



and it is the saddle point spectral density for 
/ Y[dX, A(A)2exp-7V(i-f-^)^A,2. (50) 

i 

This is 

9 , 

(51) 



Pt 



(A) = —y/r^-X^ 



with r'^{t) = AtiT - t)/T. One can compute Ut 
similarly and check that the resulting 

p+(A; t) = A (^(1 _ |)A - y^^^^ (52) 

is a solution of (^5|). 

In fact, this is the weak coupling solution for 
the YM2 sphere problem as well (found this way 
in 11), with T = A. The periodicity 9 9 + 2Tr 
of eigenvalue space is completely invisible until 
we reach r{A/2) ~ tt. This is the position space 
version of what we saw in section 2. 
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Once we reach A = n, the non-Hnearity of (|45| ) 
means this is no longer a solution. In the phase 
space formalism this is due to the constraint p < 
Pmax- One could do a calculation similar to that 
of the Wilson loop on [Hl|] to find the full 
strong coupling solution. 

In this language, the difference between the two 
problems is subtle - the plane problem 'knows 
about the compactness of eigenvalue space from 
the beginning,' while the sphere problem 'does 
not know until it hits the constraint.' It all 
comes from the boundary conditions: whereas the 
sphere problem can be completely taken over to 
non-compact eigenvalue space, the final boundary 
condition for the plane p{6) = 1/2-k only makes 
sense for compact eigenvalue space, and enforcing 
it affects the solution at all times. 

The compactness enters in a slightly strange 
way in the initial value description: although the 
original problem was free fermions, which a pri- 
ori would not know whether they were moving 
in compact or non-compact space, when we peri- 
odize the boundary condition ( |4q), we change the 
motion around = 0, because ( [45| ) is non-linear. 

The goal of this section was to study the large 
N transition entirely in terms of invariant observ- 
ables, since this language can generalize to I? > 2, 
and allows contact with a string reformulation. 
Clearly any conclusions we try to draw for D > 2 
will be speculative, but let us imagine we have a 
local continuum regulated gauge theory in D > 2. 
We need a cutofi^ and thus the dependence of ob- 
servables on length scale and coupling constant 
are not directly related. We will study it on M^. 

We can study the dependence on length scale 
by doing radial quantization of the theory, which 
turns scale dependence into 'time' dependence. 
Like the YM2 expectations on the plane this 
should be a classical boundary value problem, 
and the large radius boundary condition provides 
a way for the system to 'know' about the non- 
trivial topology of configuration space and avoid 
a transition as a function of length scale. 

On the other hand, the strong coupling expan- 
sion has no such mechanism - it only sees the limit 
g — > 00 in which the constraints, the non-trivial 
structure of configuration space, are invisible. 
Thus summing it will produce results which have 



no reason to satisfy any constraints. The same 
conclusion would apply to any string reformula- 
tion in which the state with zero strings is the 
extreme strong coupling state Wo{L) = N5lj- 
The true W{L) is produced by summing string 
diagrams, and nowhere in this formalism are the 
constraints enforced. 

YM2 on the plane is a special case because 
the observables depend only on g^A and thus the 
first argument eliminates the transition in g. |^ 
At present, all we can infer from this is that the 
second argument does not prove there will be a 
transition. 

What about a system in finite volume? Let us 
start by describing another way to understand the 
transition for S*^, developed by D. Gross and A. 
Matytsin |2^] : it is driven by non-trivial classical 
solutions in the path integral. One can evaluate 
the path integral for Z{S^) as a sum over classi- 
cal solutions, using a non-abelian generalization 
of the Duistermaat-Hcckman formula. |5^] These 
are unstable solutions contained in U{1)^ and 
constructed from the U{1) instanton. (There is 
no SU{N) instanton in _D = 2 as bundles over S*^ 
with gauge group G are classified by tti{G).) 

The U{1) instanton is the configuration Dirac 
proposed as a monopole in D = 4, restricted to 
constant r and t. As usual for instantons in large 
N gauge theory it has 0{N) action, from the ex- 
plicit N in e~^^ . However it can survive the large 
N limit because of 'entropy' - each U{1) factor 
has an independent instanton number - and in 
the strong coupling phase it does. This could 
be thought of as a covariant description of the 
'windings of eigenvalues around S^' which would 
appear if we made the classical fermion picture 
discussed here completely explicit. 

Gross and Matytsin propose that the D — A 
transition could also be driven by topologically 
non-trivial configurations, the standard instan- 
tons. In |2^, they make a one-loop estimate of 

^ In other words, a reformulation in which this result for 
the limit <; ^ 00 is reproduced in the 'obvious' way - by 
the world-sheet Boltzmann weights going to zero, as op- 
posed to cancellations between large numbers of diagrams. 
^ Conceivably, we could derive a string in D > 2 from a 
continuum formulation with no cutoff dependence or di- 
rectly from the renormalized theory, in which case this 
would also apply. 
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their effect as a function of the volume of the 
system (which controls the maximum of the run- 
ning coupling) and find a result consistent with 
no phase transition. (See also [^ ). 

From the reasoning we gave earlier, we would 
conclude instead that we generally expect tran- 
sitions as a function of volume in D > 2. At 
some critical volume, a mode will be strongly 
enough coupled to explore its full configuration 
space (e.g. < 6 < 2tt), new constraints will 
come into play, and new solutions of the collective 
field theory can exist. Now Gross and Matytsin 
are careful to state that they are arguing against 
a large N instanton induced phase transition so 
there is no direct contradiction between these two 
statements: the instanton is certainly not the 
only field configuration which can explore the full 
configuration space. 

In any case these questions must be considered 
unsettled, and the main point of my discussion is 
to outline considerations which we might be able 
to make precise given more control over D > 2 
collective field theory (or other loop equations). 
The main problem for now is to find a continuum 
regulated theory we can work with in £> > 2. 

So far I have essentially been identifying 'gauge 
string' with 'continuum strong coupling expan- 
sion.' Now this is really an assumption and if it 
turns out that the strong coupling expansion is 
invalidated by a large N or other phase transi- 
tion, we would certainly want to change it. 

What went into this assumption? The main 
element was the assumption that the state (or 
configuration - the discussion is really parallel in 
the canonical and covariant formulations) with no 
strings is the loop functional Wo{L) = N6l,i- 

As we saw, 'strings' are a reformulation of per- 
turbation theory, not necessarily in the coupling, 
but in the sense that they describe 0(1) varia- 
tions of an 0{N) initial configuration. Thus the 
choice of initial configuration to expand around is 
crucial. There is no way to describe constraints 
such as pl(0) > if we are restricted to state- 
ments expressable in perturbation theory around 
the extreme strong coupling state. This is a good 
motivation to change our zeroth order state. 

One candidate which comes to mind is the ex- 
treme weak coupling state, Wi{L) — N for all 



loops. Thinking about this leads us to ask the 
question, can standard weak coupling perturba- 
tion theory be reformulated as a theory of con- 
tinuum world sheets, if we restrict ourselves to 
computing gauge-invariant observables? This is 
an interesting question which we will not discuss 
here. In any case this perturbation theory is seri- 
ously flawed by infrared divergences coming from 
the massless gluons. These do not cancel in all 
gauge-invariant quantities and no consensus has 
emerged as to whether it can be used even in prin- 
ciple for low energy physics.^ 

Both starting points are handicapped by their 
serious qualitative differences with the true vac- 
uum of the theory. Now in principle, we could 
imagine taking any loop functional Wg{L) as our 
starting point for perturbation theory. It may 
not even be necessary for it to be the solution of 
some solvable limit of the theory. Let us write 
the true solution as W{L) = Wg{L) + 6W{L), 
and imagine deriving an expansion by plugging 
this into (^). In general we will get terms of or- 
der {dW)\ (SW)^ and if Wg{L) does not solve 
the equation, {6W)'^. We will want to interpret 
the first and second terms as leading to a 'string 
propagator' and 'string vertex' - and the third 
term as a 'string source,' which might be mod- 
eled as holes in the world sheet, but much better 
as additional pieces of world-sheet with an action 
chosen to produce the correct dependence on L. 

We could also change the association of small 
fluctuations of the loop functional W{L) with 
'strings,' as long as we preserve the locality of 
the resulting world-sheet theory. 

To get new starting points, we need to have 
some non-perturbative way to construct loop 
functionals. Now the space of loop functionals is 
sometimes said to have little structure: one just 
specifies W{L) for each L satisfying |VF(i)| < 
N VL. This is not true and in fact the space is 
not easy to describe in D > 1. What is easier to 



*What is better understood is to introduce an IR cutoff 
by hand, either by Higgsing, or by putting the system in 
a small finite volume in which the coupling does not get 
strong. The classic st udy of the large N weak coupling 
perturbation theory is [ETf . More recently Yang-Mills the- 
ory (at finite A'^) has been shown to be rigorously definable 
in finite volume. |3ql 
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describe are master fields. 

4. Master Fields 

Although we argued that taking advantage of 
factorization required us to restrict attention to 
U{N) invariant quantities, in YM2 it was very 
useful to talk about the holonomy for a given loop 
and its spectral density. 

The higher-dimensional generalization of this 
is the 'master field,' a single gauge configu- 
ration A^{x) defined to reproduce all Wilson 
loop expectation values by evaluation: W{L) — 
tr Pexpi J^Afj {x)dx^. Since the only assump- 
tion is factorization, the idea is also valid for large 
N vector and matrix models, and the master field 
has been worked out for some vector models. | |3^ 

Physicists have not had much success in writing 
explicit master fields for higher dimensional gauge 
theory or matrix models. It turns out, however, 
that recent mathematical work, specifically in the 
field of operator algebras, allows this to be done. 

The most extensive work has been done by D. 
Voiculescu and collaborators. I will only give an 
introduction here, leaving out many interesting 
results, which are described very clearly in |52|. 
More recently, I. Singer has shown how to con- 
struct the master field for YM2. Q 

The problem we discuss here is to compute 
leading order expectation values under integrals 
over several hermitian matrices Mi = . Let 
the index range over the set S. Using factoriza- 
tion, the general case is 

- ^ / n ^"''^^^ ^'"""^ j^^'- ^ (53) 

where W is an arbitrary 'word' or product of ma- 
trices 

W = M21M2^^M2I ... (54) 

We will restrict ourselves to the case 
y=^try„(M„), (55) 

in other words no coupling between the matrices. 
Thus expectation values (tr M^) are determined 



by the usual saddle point analysis and spectral 
densities pa{Xa), and this is hardly a 'c > 1' ma- 
trix model. Nevertheless we have retained an as- 
pect of the c > 1 problem: the total number of 
words of a given length I, \Wi\, grows exponen- 
tially with I, and it is non-trivial to compute their 
expectations just given the Pa{Xa). 

Still, it is not very hard to do it ~ one can con- 
sider words with successively larger m in (^), 
and write a chain of Schwinger-Dyson equations 
determining the case m in terms of lower order ex- 
pectations. Nevertheless this is rather unwieldy 
and one can ask for a simpler description of the 
result. This could be provided by a 'master field' 
- a set of single matrices "Ma = (Ma)" which 
we could regard as the saddle point for the inte- 
gral (|3|) . Now for each of the individual matri- 
ces, there is certainly a diagonal matrix with this 
property, but we cannot simply use this set: e.g. 
consider V ^ tr + tr for which 

(tr A^B^) =N (56) 

while at leading 0{N) 

(tr ABAB) = 0. (57) 

Now an integral like (^3|) is not literally domi- 
nated by a saddle point for the Ma. If we write 
Ma = UaDaU^, clcarly fluctuations of the rela- 
tive angular degrees of freedom UaU^ are com- 
pletely unsuppressed. On the other hand, the 
graphical argument for factorization can be made 
rigorous, and it is not clear why the heuristic ar- 
guments for the master field should fail. 

The reconciliation of these two statements is 
that we can contruct a master field which behaves 
as if it were the saddle point. Let us do it for the 
case Va{Ma) = ^tr M^. 

We define the 'free Fock space on S'^ to be a 
Hilbert space with orthonormal basis vectors la- 
belled by an ordered list of zero or more elements 
aj € S. The elementary operators acting on this 
are la and Z* defined as follows: 

la\an,an-i, ■ • ■ , ^i) = \a, an, ■ • ■ , ai) 

^a \ 1 ^n~l i ■ ■ ■ i ^a.an l^n— 17 ■ ■ ■ : ^l) 

^In this section we will always use the word 'free' in its 
mathematical sense, to mean a group or algebra with no 
non-trivial relations, rather than its physical sense. 
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L\) = 0. (58) 

They are similar to bosonic creation and annihi- 
lation operators but with two differences. First, 
they are free - the product of two operators as- 
sociated with a and (3 ^ a satisfies no relations. 
Second, the usual symmetry factor for bosonic 
harmonic oscillators is absent. Thus, we have not 
[/, r] = 1 but instead 

E[^-^^] = l) (I- (59) 

On this Hilbert space, we define an operator 
Ada associated with each M^, and a 'trace' tr . 
For the gaussian integral V = 

Ma = L + ll- (60) 

The 'trace' is not the standard one (which does 
not make sense here) but is defined by 

tiA^{\A\). (61) 

It is a trace in that it satisfies the axioms for a 
trace, for example tr [A, B] = 0. (This is not yet 
obvious.) 

One can see graphically that the construction 
works, 

ti Ma, . . . Ma„^ = (tr Ma, . . . Ma J, (62) 

by expanding the product of {lat + ^a ) 2™ 
terms and associating a planar diagram with each 
non-zero term. Consider the example in figure 4. 
We start with Ma^=a\) which makes the state 
I a), and we interpret this as adding a line with 
the label a to the diagram. With Ma^_i, there 
are in general two possibilities: we can always 
act with la^ -, to create a new line labeled am-i, 
and if am-i = Oim we can annihilate the exist- 
ing line with la^-i as well. The 'free' nature of 
the Fock space precisely reproduces the planarity 
constraint on the diagrams. After applying Mq,,, 
we must be left with no lines to have a matrix 
element with ( | . This proves ( |6^ ) , and since tr is 
a trace, so is tr . 

The diagrammatic argument is a good demon- 
stration but to gain insight from the use of master 
fields we need new concepts. The central concept 
in Voiculescu's work is that of 'free probability 




Figure 4. Each non-zero term in the expansion 
of (tr MfcM^MfjM^) corresponds to a planar dia- 
gram. 



distribution.' This is analogous to the idea of 
independence in probability theory: a set {xa} 
of commuting random variables is independent 
if the joint expectations 4>{fi(xi) . . . fn{xn)) = 

(l>l{fl{xi)) . . . (j)n{fn{Xn))- 

Let us regard our matrix model expectation 
values (tr W) as giving the joint expectations 
(/)(VF) for words constructed from a set {Ma} of 
non-commuting random variables. We normalize 
0(1) — 1. The distribution is free if 

0(/,(M„J) = O i = 1,2, . . .m and (63) 

ai 7^ aj+i Vi 

imply 

0( fl{Ma,) ... fm{MaJ ) = 0. (64) 

Just as for independence, given the individual dis- 
tributions this assumption completely determines 
the joint distribution. The general expectation 
0( fi{Mai) ■ ■ ■ fm{Ma^) ) cau bc computed by 
induction. Let 5.(M„J = /,(M„J - (j){f,{Ma^)) 
so that 4){gi{Mai)) = 0; then by substituting 
fi = 4>{fi) + gi for all i, distributing and using 
( |64| ) on (j){gi . . . gm), we can reduce it to a sum of 
terms with fewer components. 

Expectations of the large N limit of decou- 
pled matrix integrals with ( p5| ) are free. This can 
be proven diagrammatically, by induction on the 
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number of components: given (63), to get a non- 
zero expectation we must take diagrams with a 
hne connecting two components, and this sphts 
the diagram into a product of expectations with 
fewer components. It has been proven in for 
a broader class of decoupled measures. 

Freeness can be implemented directly in the 
master field: given the individual master fields 
for the individual matrix integrals, which we al- 
ready know, the joint master field is their 'free 
product.' This can be defined algebraically us- 
ing the ideas we already described: the essential 
point is that Mi and Mj satisfy no relations. 

To fit this into a larger context, we will change 
our language slightly: we should think of the mas- 
ter field as a representation of the abstract alge- 
bra generated by the M^'s, and ^(^) as a linear 
functional on this algebra. Many of the usual 
ideas of representation theory, such as unitary 
equivalence, are directly relevant here. As we will 
see later, this language is also well suited to the 
gauge theory case. 

The relevant algebra for a hermitian matrix in- 
tegral is the algebra of functions of a single real 
variable A. For mathematical purposes of course 
we would want to be more precise. I will not 
try to be so precise, but at least for our cultural 
benefit let us look at the algebras which come 
up. For a one-matrix distribution, M is repre- 
sented by the operator of multiplication by A (its 
spectral parameter), and (t>{f ) = / dA p(A)/(A). 
This only sees a subset of M, the spectrum of 
M (support of p), so we could equally well use 
the algebra of functions on SpecM. One natural 
object to consider is the algebra of bounded con- 
tinuous functions on SpecAf. This can be made 
into a C*-algebra by using the sup norm. (For the 
defining axioms of a C* -algebra, see for example 
[Q.) C*-algebras are the most tractable general 
class of algebras, and this is quite useful. It is a 
theorem that every commutative C*-algebra is an 
algebra of bounded continuous functions on some 
space M , and one can recover the topology of M 
from the algebra C^{M). In matrix models, it can 
be relevant to the physics to know the number of 
components of SpecM (e.g. one-cut v.s. two- 
cut matrix models) and their topology, and per- 
haps analogous concepts of the 'non-commutative 



topology' of the algebras arising in multi-matrix 
models or matrix field theories will eventually be 
relevant. 

Since we have a measure, another natu- 
ral algebra to consider is that of the L°° 
fimctions, bounded in the norm ||/||oo = 
limp_^oo {Jd\ p{X)\f{\)\P Y'^. This is both a 
C*-algebra and a von Neumann algebra. These 
functions need not be continuous and this algebra 
retains very little information about the topology 
of the underlying space. 

Usually for us SpecM is a bounded set /, so 
we can restrict attention to C°(/) or L°°{I). 
This last is isomorphic as an abstract algebra to 
L°°(S'^), and by Fourier transform to L(Zi), the 
group algebra of Z. | Now wc can appeal to a 
result of which I hope will sound intuitive: 
the free product of group algebras is the algebra 
of the free product of groups. The free product 
Z * Z . . . {n times) is the free group on n gener- 
ators, Fn- Thus the von Neumann algebra which 
arises in an n-niatrix problem is L{Fn). 

It is conceivable that other algebras could come 
up. For example, in the large N limit of Chern- 
Simons theory, ||l^ one expects Wilson loop ex- 
pectations on S X / (for S a Riemann surface) 
to give a functional on L(7ri(S)). In general, we 
know that we will get a 'type IIi algebra in the 
Murray- von Neumann classification.' This just 
means that there is a trace taking values (on pos- 
itive operators of norm one) in [0,1]. One also 
has a notion of 'factor,' an algebra with trivial 
center, which implies that this trace is unique. 
This may sound attractive for our purposes but 
at this point one starts getting into subtleties, in 
particular involving the choice of norm used in 
defining the algebra. There are many definitions 



®The group algebra of G, C(G), is the algebra of linear 
combinations of group elements ^2g£G ""^^^S with multi- 
plication induced from the group law, and g* = ■ It 
is essentially the convolution algebra of functions on G. 
Perhaps the simplest precise version of this is to take all 
functions on G (with an invariant measure). For 'nice' 
groups such as compact Lie groups and abelian groups, 
the representation theory can be used to turn convolution 
into multiplication on a dual space, as we were doing in 
section 2, and for the group Z, this is just Fourier trans- 
form - it is a good exercise to check this if you haven't 
seen it. 
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of C(F„) (for n > 2): some are factors (e.g. the 
'reduced' algebra), and others are not (e.g. the 
'full' or 'universal' algebra). The only point I will 
take from this is the obvious one that one needs 
to be careful to check that one's definitions are 
appropriate. 

Returning to a single matrix, another descrip- 
tion of the same representation which makes con- 
tact with our free Fock space is to write a series 



L + 



(65) 



satisfying tr — J dXp{X)X''. Clearly such a 
thing exists in general, because the fc'th moment 
is determined by c„ with n < fc— 1. The free prod- 
uct of these representations can then be taken by 
letting these operators act on the free Fock space. 

One can show that the c„'s are determined 
as follows: given the resolvent F{z) = (t){ ^\.j ) = 

J '^z-x^ 1 compute the inverse of F{z) under com- 
position, K{F{z)) = z. Then 



K(x) = -+ nix) = i + V c„a;" 

X X ^ — ' 

n>0 



(66) 



For example, the semicircle distribution has 
F{z) = ^(z - Vz^ - r2) giving n{x) = ''^x. 

Thus, the operators ( |65| ) acting on free Fock 
space are the master field for any decoupled set 
of matrix integrals. What can we do with this? 
A problem treated in is the following: find 
the moments 0((X]a -^a)")' O'' equivalently the 
spectral density for Ma. In ordinary proba- 
bility theory, the resulting distribution would be 
the convolution of the independent distributions; 
thus the new concept can be called 'free convolu- 
tion' (or 'additive free convolution'). 

It is clear that if we add two free variables, the 
result will be free with respect to the rest, and the 
addition is commutative and associative. Thus we 
need only consider binary free convolution. The 
result which determines it is the following: the 
operators 

have the same distribution as 



This is not hard to see, by explicitly writing out 
a term like (|(Mi + M2f\) = + h + ■ ■ 
and noting that any place an annihilation opera- 
tor appears, exactly one of li or I2 will contribute. 
Thus 



7^l+2(a;) =7^l(a;) +7^2(a;) 



(69) 



For example, the free convolution of semicircles 
with radius is a semicircle of radius = X) ■ 

In an analogy is developed between the 
Gaussian in ordinary probability theory and the 
semicircle law in free probability theory. For ex- 
ample, there is an analog of the central limit 
theorem: given a free family of distributions (pi, 
with (j){Mi) = 0, all higher moments bounded, 
and lim„^oo ^ T,i=i HMf) = > 0, the sum 
Si converges on a semicircle with radius 
r. Such results might help explain the ubiquity of 
the GOE and GUE in statistics of energy levels 
of random Hamiltonians. 

So far we have been considering free products of 
one matrix distributions but the concept is more 
general: for example, we could take the free prod- 
uct of (j) with arbitrary dependence on Mi and 
M2, with -0 depending on M3. In full generality 
we could talk about a free family of subalgebras 
Ai of the full algebra of ordered products: in (|64[), 
the successive components would be required to 
belong to distinct subalgebras. 

Let us use these techniques to give another 
derivation of (^|) , here for quantum mechanics of 
a hermitian (rather than unitary) matrix. (This 
is very similar to |Q ) . We start with a discretized 
path integral: 



Z = Jf[dM{^) e-^S.'^ W^+i) 



(70) 



with t = eL. Given an initial condition M{0), 
this determines 



Af(i) =A/(0) + ef]rK*) 



(71) 



Mi+2 = l + }_^ici. 



C2.n)l* 



(68) 



where the 77(1) = M{i) — M{i — 1) are free ran- 
dom variables with semicircular distribution. Let 
Pt{X) be the resulting spectral density at time t, 
and define Ft and Kt from it as above. We can 
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compute Kt using ( |69| ) in terms of the initial con- 
dition Kq{u): 

Kt{u) ^ Kq{u) + eLu (72) 

Changing variables from u to ^ = Kq(u), 

Kt{Fo{0)=^ + tFo{0 (73) 

Let us define z{^,t) = Kt{Fo {£,)). Applying 
Ft{z) to both sides of this definition, we have 



Ft{z) = FoiO 



(74) 



and we have reproduced (^ ) with Ft{z) = 
p+{z,t). Finally, substituting into ( [73| ) produces 
(^), the general solution of (^. 

Notice that this derivation only required an ini- 
tial condition and produced the Hopf equation di- 
rectly for the resolvent tr {z — AI)^^. Perhaps we 
have found a better explanation for the condition 
P± ~ 1/^ we found in YM2 on the plane? That 
result was almost the same but with periodic 6 - 
can we derive it in this spirit? 

One can imagine calculating the master field 
for the holonomy in figure 5 by summing graphs 
in axial gauge. This would produce an evolution 
similar to (|7l|), but operating by multiplication: 



[/(< + £) = e'^(*'0)C/(i)e-*^(*'^) 



(75) 



with the {A{t, x)} at different times free with re- 
spect to each other. This should be calculable us- 
ing 'multiplicative free convolution,' another con- 
cept developed in |52[ . 

A direct field theory application of this is 
to Gaussian matrix models (free in the physics 
sense). The Z3-dimensional matrix field theory 
with (Euclidean) action 
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(76) 



is a free product of semicircular distributions for 
each mode A{k) + iB{k) = J d^a;e*'=^M(x). 

It is not clear whether freeness plays any direct 
role in non-trivial higher dimensional field theo- 
ries, but it may be possible to use free distribu- 
tions indirectly to construct their master fields. A 
construction of this type was given by J. Green- 
site and M. Halpern: let us discuss it briefiy. |22] 
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Figure 5. Time evolution in YM2. 



A general field theory can be formulated using 
stochastic quantization: |Tl| ] the t 00 limit of 
the solution of the Langevin equation 



d_ 
d'T 



0(x) 



dS 



d(t){x,T] 



+ 7]{x,t), 



(77) 



where 7y(x, r) is a Gaussian random variable 
with correlation (77(0:1, Ti)?7(a;2, T2)) — 2S{xi — 
X2)S{ti — T2), reproduces quantum expectations, 
{4'{x)...)rj = {(j){x)...)^. 

For (j) a hermitian matrix, the ri[x,T) be- 
come free random variables in the large N limit. 
(Greensite and Halpern represented them using 
explicit Gaussian matrix integrals.) Thus we can 
regard ( [77| ) as a 'free non-linear PDF' for which 
the r — > 00 limit of a solution is a master field. 
This is an equation with some similarity to (|7l|), 
and perhaps one could find free analogs of the 
existing techniques to solve non-linear PDE's (a 
particularly interesting idea for integrable theo- 
ries). 

Let us at least solve the linear case: for (|7^) , we 
can Fourier decompose cf) = M , and proceeding as 
we did with (Fll) for each mode, find that 



d 

^Kr{k- u) ^ -{P + m^)Kr{k; u) 

OT 



u. 



(78) 



The solution is Kr{k;u) = u/{k^ + m^) + 
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Cue~^^^'^"^^^'^ which for r ^ oo corresponds to 
the semicircle we mentioned earher. 

Another possible use of free expectations would 
be as ansatzes or as starting points for per- 
turbation theory. An example of something 
we can get this way is a gauge theory loop 
functional which interpolates smoothly between 
'weak coupling' behavior of the spectral density 
at high momenta and 'strong coupling' behav- 
ior at low momenta. It is just defined by the 
integral / dAexp - / d^fctr /(fc, Af,{k)A^'{-k)) 
(with A^{—k) = A^[k)*). f{k,z) can be chosen 
arbitrarily and to get the required property we 
would take it to interpolate between /(fc, z) k'^z 
at large k, and (say) f{k, z) ~ fc^ cos ^/z at small 
k. 

The ansatz is rather ugly, especially with its 
ad hoc gauge fixing. Nevertheless it does define a 
master field and all Wilson loop expectations can 
be computed from it. A possible apphcation for 
it would be as a 'building block' for more general 
loop functionals, to which we turn. 

There is a more general construction, by which 
if one is given all expectation values of invari- 
ant operators, one can construct a master field 
which reproduces them. This has been done by I. 
Singer for a continuum gauge theory, specifically 
YM2. Let us discuss an analogous construc- 
tion in lattice gauge theory. (One can also adapt 
it for matrix models.) In this case, the gauge 
holonomy is a representation of the groupoid of 
paths on the lattice.]] For the present discussion 
(this can be generalized) we fix an 'axial' gauge 
by choosing a spanning tree T of the lattice, and 
setting these links to 1. The holonomy is then 
determined by its values on the graph formed by 
contracting T to a point t; this graph is a 'bou- 
quet' of circles attached at t and paths starting 
and ending at t generate a free group algebra. 
Thus W{L) is a functional on the free group, and 
we define its values on a general element of the 
algebra by linearity. 

The construction of the master field is the GNS 



'^A groupoid is like a group, but the multiplication law 
need not be defined between all pairs of elements. |^ Here 
the multiplication pi o p2 is defined if the end of pi is the 
start of p2- In the associated algebra, a multiplication 
undefined in the groupoid produces zero. 



(Gel'fand-Naimark-Segal) construction and 
once one has phrased the expectation values as a 
linear functional on the group algebra, it applies 
straightforwardly. The regular representation tt 
of the algebra A acts on a vector space H = Ahy 

TT{a)\b) = \ab). (79) 

Essentially, we want to make this a Hilbert space 
with the inner product 

{a\b) = W{b*a) (80) 
and define 

tr 7r(a) = (/|7r(a)|/). (81) 

(We are labelling vectors by the associated ele- 
ment of A, and / is the identity element in the 
group.) 

The master field will be the representation 
7r(L). Its defining property is that tr 7r(L) = 
W{L); this is almost a tautology and in this sense 
we only 'repackaged' the information of the loop 
functional. But there is more to check, because 
we need the inner product to be positive definite, 
TT to be a unitary representation, and tr to have 
cyclic symmetry. 

The first condition clearly will not hold for just 
any W{L): it must satisfy a positivity condition. 
Now positive definiteness of the inner product 
would follow from W{a*a) > Va, but it is clear 
that semi-positive definiteness 

W{a*a) > Va (82) 

is all we can expect in general. Let's reconsider 
the D = 1 case: loops are classified by winding 
number and their group is Z; the positivity con- 
dition is 

^W^„~™/„/™>0 V/„. (83) 

We can turn the convolutions into multiplication 
by considering p{9) — ^^Wn^™^ and in 
terms of which ( |83| ) becomes 

1 ddp{Q)\!m^ > (84) 

which is our old positivity condition. 
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One can define a positive element p of an oper- 
ator algebra as one whose eigenvalues are positive 
in any representation, or more abstractly by re- 
quiring its spectrum (this is defined without ap- 
peal to any representation as {\\{p — A)^^ ^ A}) 
to be contained in the positive real axis. An ele- 
ment p = a*a will be positive, and it is a theorem 
for C* algebras that any positive element can be 
written in this way. One refers to as a positive 
linear functional if it is positive on the positive 
elements. 

We now deal with the semi-definite case by 
finding the ideal of elements {2;|VF(z*z) = 0} 
and quotienting the space Ti by this ideal. In 
the one matrix case, the quotient is the algebra 
of functions on Supp p. The inner product ( |80| ) 
is positive definite on the quotient and there is 
a unique completion of the quotient which is a 
Hilbert space. The result is a unitary representa- 
tion. Unitarity u* — follows because it was 
true in the abstract algebra. 

Finally, tr [a,b] = (l|[a,6]|l) = W{[a,b]) = 0. 
Thus the representation tt{L) is a master field. 

Not every loop functional has a master field, 
because of the positivity requirement. A (uni- 
tary) quantum field theory vacuum, however, 
will. First, the contribution from any point in 
field space is positive: by definition, a positive 
L must have positive eigenvalues in any repre- 
sentation, and thus tr U{L) (for positive L) will 
be positive for any U. Then, expectation values 
under a functional integral with positive measure 
will be positive, so these will be associated with 
master fields. 

If we need to deal with a non-positive loop func- 
tional (say a perturbation of the ground state) , we 
can construct a 'virtual master field' which is a 
pair of fields Ui, U2 such that W{L) = tr Ui{L) - 
tr U2{L). It can be proven that this exists for any 
real W{L)^ and that a similar construction with 
four Ui can reproduce any complex functional. 

If we regard the large N theory as being defined 
by a loop equation, e.g. the Migdal-Makeenko 
equation or equations derived from a collective 
Hamiltonian, we should regard the positivity con- 
straint as an additional element of the definition. 
It is non-trivial, non-perturbative structure of the 
configuration space. To some extent it is indepen- 



dent of the loop equation, and when we solve the 
equation we may need to enforce the constraints 
by hand, as we discussed in the previous section]^ 
These inequality constraints are rather simpler 
than the Mandelstam constraints on loop func- 
tional constructed from gauge fields at finite N 
~ for example, there are cases (loop functionals 
near the strong coupling limit) in which they are 
irrelevant. It was known in earlier work that such 
constraints exist and important, and the de- 

scription given there had some similarities: when 
the change to loop variables is non-singular, the 
metric on configuration space {hij of section 3) in 
these variables, 

17(L,L') = 5]tr [ElW{L)][ElW{L')] (85) 

L" 

will be positive definite. The constraint (|2|) is a 
bit stronger than ( |85| ) (in D > 2) because it does 
not have the sum over L" . It is simpler to think 
about, most importantly because it is obviously 
linear. It is also clear that positivity requires pos- 
itivity for any subgroup of the group of loops. 

Let us illustrate it on a 'figure eight,' a bouquet 
of two loops li and I2 with holonomies Ui and U2- 
Besides a list of the W{L), an equivalent form of 
the data is to specify 

W{F,Hf ...F-ll-)= (86) 

The simplest new case of (^) is to let a = 
Srn n Cm,nT'2 • thcsc Can be Combined to show 
that /9i(^i,0i) > 0. We could also de- 
rive this from a spectral decomposition Ui = 
/ de p^{e)P^{e) over projections P^{ef = P^{0), 
using tr P P' > 0. However the pattern does 
not generalize: the next constraints, derived from 
" = Y^i.m,n ci,m,nli^Tli^ ^^e too Weak to forcc P2 
to be positive. It is not hard to construct exam- 
ples in which it is not, using projections such that 
tr P P' P" < 0. 

Although this does illustrate that (|2|) contains 
new information, this particular explicit form 

* To complete the discussion, we should prove that every 
positive loop functional can be produced by some path 
integral and thus there are no further constraints. We 
have only done this for small modifications of the strong 
coupling vacuum and this is an interesting open problem. 
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does not look too useful. Another illustration of 
the difficulties of working explicitly with (^2|) is 
the following: there are functionals on a bouquet 
of k loops which satisfy the constraints defined 
in terms of Z < fc loops, but violate constraints 
involving k loops. 

A loop functional constructed from a master 
field automatically satisfies the constraints, and 
(as was noted in previous work) the simplest de- 
scription of the configuration space is the space of 
all master fields up to gauge equivalence, or (what 
is the same thing) the space of unitary represen- 
tations of the free group algebra up to unitary 
equivalence. 

An important issue in working with D > 2 
large N theories is the choice of representation 
of loop functionals - for both analytical and nu- 
merical purposes we need truncations which al- 
low approximating a loop functional with a finite 
amount of information. If constraints are satu- 
rated, it is clear (already in _D = 2) that specify- 
ing the functional on all loops up to a given length 
is inappropriate. One alternate proposal in ear- 
lier work was to work with a finite-dimensional 
approximation to the master field and derive the 
loop functional from it. Another style of represen- 
tation is to build loop functionals up from simple 
functionals. Since the space of loop functionals 
really is a linear space, it makes sense to propose a 
basis for loop functionals. One might try to find a 
basis of 'extreme points' in the space of function- 
als, from which any positive functional could be 
made by superposition with positive coefficients. 

Simple operations which preserve positivity are 
addition of loop functionals, and loopwise multi- 
plication (Wi X W2){L) = Wi{L)W2{L). For- 
mally, the corresponding master fields are just 
the direct sum and direct product of those for Wi 
and W2. Although we have not precisely defined 
continuum master fields, it is interesting to note 
that in the continuum, if Ai and A2 are solutions 
of the classical Yang-Mills equations, both direct 
sum and direct product will produce a new solu- 
tion. This explains the observation of A. Migdal 
that the left hand side of the Migdal-Makeenko 
equation (^) is a first order linear operator on 
loop functionals. po| 

There is a similar idea of 'phase space master 



field' which reproduces expectation values of in- 
variant operators made from A{x) and E{x), or 
for matrix models M[x) and its conjugate T\{x), 
as used in . Issues of operator ordering make 
this rather more subtle to study, but we intend 
to discuss it in future work. 

Let us discuss a simpler question raised in the 
last section: what is the higher dimensional ana- 
log of the Virasoro algebra of matrix quantum 
mechanics, and the corresponding group? The 
generators are known: 

n{!,x) = iv f{M)n{x) (87) 

where f{M) can be an arbitrary word in the M's 
(but not the n's).| Their algebra is easy to com- 
pute as is their action on the invariant observables 
c^{W)=tr M{xi)...M{xmy- 

[Uif,x),^iW)]=Y,S-,-. (88) 

i 

tr M(xi)...M(x,_i) / M{x^+l)...M{x„^) 

Can we generalize to D > 1 the identification 
of this as the infinitesimal diffcomorphism 6X ~ 
/(A) ? In terms of the master field, the general- 
ization is 

6 Mix) - f{M) (89) 

and we should think of this transformation as an 
'infinitesimal free diffcomorphism.' 

It is an infinitesimal automorphism of the free 
group algebra C{Fn): it is linear and a derivation, 
and SM+ = {SM)+ implies it is an automorphism 
of the C* algebra. Thus the corresponding group 
is the automorphism group of the free group alge- 
bra (or possibly the connected component of the 
identity) , and the coherence group is a semidirect 
product of this with the additive group in the al- 
gebra. 

Returning to configuration space, a simple but 
important master field for gauge theory repro- 
duces the extreme strong coupling limit, the loop 
functional W{L) = Slj. It is just the regular 
representation with the conventional inner prod- 
uct {Li\Lj) = Si J in a basis of loops. It saturates 

^ We can assume f{M) = f{M)'^ without loss of gen- 
erality. Strictly speaking, we should define !!(/, x) = 
^tr (/n -|- n/), a self-adjoint operator, but this is not im- 
portant for the present discussion. 
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no constraints: by adding sources to the func- 
tional integral we can make small variations of 
the W{L) for every L independently. 

The first observation to make is that this mas- 
ter field is a representation on a very large Hilbert 
space, not just infinite but 'more infinite' than 
even a quantum field theory Hilbert space. What 
I mean by this is that in terms of the natural 
grading, namely by word length, the degener- 
acy of states grows exponentially as {2D — 1)' . 
(For quantum field theory in dimension D and 
with finitely many fields, the degeneracy of states 
grows as exp E^-^^^.) It is the same linear space 
in which variations of the loop functional live 
and we can (by analogy with string field the- 
ory) also think of it as the 'first quantized string 
Hilbert space.' From this point of view, expo- 
nential growth is the expected result for a string 
theory. However there is an important difference 
with the better understood case of fundamental 
string theory: although there we also have expo- 
nential density of states in terms of space-time en- 
ergy, the world-sheet grading is not by space-time 
energy but by world-sheet 'energy' n = a'rn^ , and 
thus world-sheet densities of states have the fa- 
miliar exp y/n behavior. 

Now the grading which is relevant in string 
perturbation theory is that given by the world- 
sheet Hamiltonian, and we don't know this for 
gauge strings. (We do know it in light-front quan- 
tization, at least if we call the states tr Oi\0) 
'strings,' and there the world-sheet Hamiltonian 
determines m^.) Thus it is not yet clear if this has 
relevance for physics. It turns out to be very rele- 
vant for mathematics, however. Consider the fol- 
lowing question: since we found that the extreme 
strong coupling master field is simply the regu- 
lar representation of the group algebra, why not 
follow the normal procedures of group represen- 
tation theory and decompose this representation, 
which for compact Lie groups contains copies of 
every irreducible representation? 

In fact this representation is irreducible. Many 
standard statements of representation theory 
break down for such large groups. Technically, 
free groups are not amenable, a statement whose 
definition I leave to |2^,^. Apparently, little is 
known about the space of all of its representa- 



tions. 

5. Parting thoughts 

A lot of work on large in D > 2 has con- 
centrated on trying to solve models or at least 
develop reliable numerical procedures for treat- 
ing loop equations. This is an attractive goal, 
not just for gauge theory but in my opinion even 
more so for matrix models and string theory, be- 
cause there the main difficulty we face - the large 
number of degrees of freedom - is not an artifact 
of our approximation but an essential feature of 
the problem we are studying. 

We appear to be far from this goal, and I 
believe that a certain amount of mathematical 
ground work will be required first, because the ob- 
jects we are working with are so unfamiliar. Thus 
I want to conclude by proposing a simpler appli- 
cation of the loop space formalisms than 'solving' 
a large A^ field theory: to be the starting point for 
new perturbative expansions, and to better un- 
derstand the expansions we have. These are prob- 
lems which require some ability to work with the 
theory non-perturbatively, may well lead to prac- 
tical calculational techniques, and could give us a 
better starting point for more ambitious work. 

The prototype for 'gauge string theory' is the 
strong coupling expansion. The recent work sug- 
gests that this old formalism can be greatly im- 
proved, most notably by formulating it in the con- 
tinuum. We need to continue it to weak coupling, 
but if it is qualitatively correct, as perhaps it 
could be, there would be no fundamental barrier 
to doing this. However, to use it one must justify 
not only the expansion around infinite coupling 
but also the assumption that the difference be- 
tween the starting loop functional Wo{L) — 6lj 
and the result W{L) is 'small' - not necessarily in 
magnitude, but in the sense that no new structure 
of the configuration space or action is required to 
get the correct result. As we saw, these issues 
are much simpler for large N gauge theories than 
for finite A^ gauge theories, but still non-trivial, 
and in the end this assumption was found to be 
highly suspect - there is new structure in the form 
of linear constraints on physically realizable loop 
functionals, which is not visible in strong cou- 
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pling, and probably leads to a large N transition 
analogous to those in the solvable models. 

We should think of 'string' as a perturbative 
construct, describing a small fluctuation of the 
fields, very analogous to particles in quantum 
field theory. Thus we should not try to fix the 
original string construction but instead define a 
new construction in which we expand around a 
different zeroth order loop functional, e.g. one 
which saturated the right constraints. Changing 
this could produce a very different string theory. 
(See for an extreme case of this.) 

Another, more conventional motivation for ex- 
panding around non-trivial configurations, is to 
get the physics right in the zeroth order. A dif- 
ferent direction would be to try to fix the obvious 
fiaws of the weak coupling expansion for gauge 
theory, such as the masslessness of the gluons, by 
starting it around a non-trivial background. 

To do any of this we need to better understand 
loop functionals and master fields, and I hope I 
have convinced the reader that there are inter- 
esting things which can be done in these direc- 
tions. I believe this would also be invaluable for 
matrix models of fundamental strings - after all, 
the real point to finding a c > 1 matrix model 
is to make non-perturbative statements using it, 
and surely the first step towards this is to have 
some non-perturbative picture of its configuration 
space. But these are problems for the future. 

It is a pleasure to acknowledge valuable discus- 
sions with M. Halpern, V. Kazakov, I. Kostov, A. 
Migdal, G. Moore, H. Neuberger, A. Polyakov, 
M. Staudacher, W. Taylor, and K. Woo, and I 
especially thank I. Singer and D. Voiculescu for 
discussing their work with mc. I also want to ex- 
press my gratitude to the organizers of the Spring 
School, and to the students for their interest and 
for their patience. 
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